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Abstract 



We study symmetry breaking in Z2 symmetric large N matrix models. In the planar 
approximation for both the symmetric double- well </> 4 model and the symmetric Penner 
model, we find there is an infinite family of broken symmetry solutions characterized by 
different sets of recursion coefficients R n and S n that all lead to identical free energies and 
eigenvalue densities. These solutions can be parameterized by an arbitrary angle 9(x), for 
each value of x = n/N < 1. In the double scaling limit, this class reduces to a smaller 
family of solutions with distinct free energies already at the torus level. For the double- well 
(f) 4 theory the double scaling string equations are parameterized by a conserved angular 
momentum parameter in the range < / < 00 and a single arbitrary U(l) phase angle. 
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1. Introduction 



We would like to highlight some unusual aspects of symmetric double-well matrix 
models [[l|-[|T(| concerning spontaneous symmetry breaking and the multiplicity of solu- 
tions at the same critical point. We consider two specific models: (i) a </> 4 theory with 
a symmetric double well and (ii) the symmetric Gaussian Penner model ||11|| - ||13|| . Both 
models have a Zi reflection symmetry and a standard two-band solution that respects 
it. However we obtain new classes of solutions that break the Z<i symmetry by relaxing 
the initial boundary conditions on the first two recursion coefficients for the othogonal 
polynomials. 



The single-band broken symmetry solution to the double-well (/> 4 theory is c=0 pure 
gravity, but at lower free energy there is an infinite class of two-band solutions where 
the tree level eigenvalue density is symmetric in the two wells. These solutions have the 
property that their free energy and eigenvalue density, in the planar limit, are invariant 
with respect to an infinite set of continuous parameters in the recursions coefficients. An 
analgous two-band class of broken symmetry solutions is found in the Penner model as 
well. In the double scaling limit for the </> 4 theory the degeneracy is lifted, except for 
a single U(l) rotation, and a one parameter family of solutions survives satisfying the 
Painleve-II equation with an extra conserved "angular momentum" parameter. At this 
stage the physical consequences of this degeneracy of solutions are not clear. Since the 
tree level eigenvalue density is the same for all these solutions it is possible that some of 
them could tunnel into each other with an instanton action that is lower order in N. If so, 
this would be of particular interest in the context of the Gaussian Penner model, where 
one of the broken symmetry solutions discussed is related to the c = 1 string at twice the 
self-dual radius. 



The organization of this paper is as follows: In section 2, we give a brief overview 
of symmetry breaking and the formalism for the orthogonal polynomial method. In sec- 
tion 3, we consider the consequences of relaxing the boundary condition on the recursion 
coefficients for the symmetric ^> 4 potential and the Gaussian Penner model, and classify 
all two-cut solutions in the planar (N — > oo) limit. In particular we exhibit the class of 
solutions that have the same tree level eigenvalue density and free energy. We discuss 
some correlators that distinguish between the various solutions and a numerical approach 
to investigating finite N solutions. Section 4 discusses the double-scaling limit of the free 
energy for both models emphasizing that one gets an expanded class of double-scaling 
solutions. The discussion and conclusions are given in section 5. 
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2. Overview of Problem 



It has been suggested that the singular behavior of the tree level eigenvalue density 
near the edge of the cuts determines the critical behavior of the matrix model in the fol- 
lowing sense |14]] || |§ : One identifies the polynomials in the matrix variable <p or scaling 
operators O n (O n 's are traces of the polynomials) which produce a particular kind of singu- 
larity (labeled by n) near the edge of the cut or cuts. Knowing the O n 's, one then considers 
Z{t) = f d(J)exp[— tnO n (4>)] and shows that Z{t) satisfies a certain hierarchy of equa- 
tions. This hierarchy depends, then, only on the kind of ensemble of matrices considered 
(hermitian, antihermitian, unitary, etc.), the class of singularities of the eigenvalue density 
allowed (e.g., single-cut density with multiple zeroes coalescing at the edge, two cuts col- 
liding and sandwiching zeroes in between, etc.), and any symmetry of the potential that 
restricts the class of eigenvalue densities considered (e.g., restriction to even perturbations 
in Z-i symmetric 2-cut models) jlQ j. 

Orthogonal Polynomials 

It is generally believed for potentials V((f>) = t n O n that are bounded from below 
that the orthogonal polynomial method [|15| uniquely fixes the solution for the free energy, 
correlators, etc. Here we will show that orthogonal polynomial method actually allows one 
to construct a whole class of closely related solutions. To clarify this further, we remind 
the reader of the precise condition for a unique solution for the free energy in terms of or- 
thogonal polynomials. Consider the partition function, Zn = J d(p e - NtrV ( ( t > )^ where (p is 
an iV x N hermitian matrix and V((p) is a real potential. The integral is expressed in terms 
of a set of orthogonal polynomials P n (x), J dx P n {x) P m (x) e~ NV ^ = h n 5 nm , normalized 

by the convention that the leading term for P n (x) is x n , P n (x) = x n -fcix^ 71-1 -' + 

(Note that this convention sets Pq(x) = 1 or ho = J dx exp[—NV(x)].) Since these orthog- 
onal polynomials can be iteratively determined by Gram-Schmidt orthogonalization, the 
exact free energy, = logZ^ = ^2 n =o ^°Shm as well as all thermodynamic averages are 
uniquely determined. 

Instead of actually finding these orthogonal polynomials, one in practice uses re- 
cursion relations for the coefficients R n and S n in the expression, xP n (x) = P n+ i(x) + 
S n P n (x) + R n P n -i(x). Once the R^s are known, the free energy can be found by using 
the fact that R n = h n /h n -\. It is convenient to introduce a self-dual orthonormal basis 
\n), where (x\n) = P n (x)/y/h~^, and \x) are eigenvectors of the operator <f> with eigenvalues 
x, satisfying the normalization (x'\x) = exp(NV(x)) 8(x' — x). Matrix elements of the 
operator <fi in this new orthonormal basis are directly related to the recursion coefficients 
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by (m\(f)\n) = ^/ R m 5 m:n+ i + S n 8 m ^ n + ^ R n b~ m>n -i. In terms of cf>, the recursion relations 
for R n and S n can be expressed in operator notation, 

n/N = y/R^(n-l\V'{4>)\n) 

= (n\V'{$)\n). ^ 

Once initial conditions are specified, a unique solution for R n and hence the normalization 
for each orthogonal polynomial, h n , can be found by iteration. For example, in the case 
of the (J) 4 model, Eq. Q is a pair of coupled two-term recursion equations for R n and 
S n which require four inital conditions: the numerical values 

S = (0\4>\0) = ho' 1 j dxxe~ NV{x) 

J , (2-2) 

R 1 = (0|<^ 2 |0> - (0|^|0) 2 = ho' 1 / dxx 2 e~ NV{x) - Si 

and the trivial values Rq = and S-i = 0, which are independent of the potential. Given 
these values, all other coefficients are given iteratively as rational function of So and R±. (A 
more elegant formulation would introduce a single sequence of coeffients Ck oc (0|(/> fc |0) + .... 
with a single four-term recursion relation Ck = F)-(Ck-i, Ck-2, Cfc-4), where Ck = 

Rk/2 f° r k even and Ck — S^k-i)/2 f° r k odd.) 

Symmetry Breaking 

This formalism should make it clear that there is no ambiguity in defining the matrix 
models for all finite N, assuming of course that the potential is bounded from below and 
that the integrals defining the inner product are finite. However typical of all statistical 
mechanical problems this does not imply that we know the correct way to take the ther- 
modynamic (or in this instance large N) limit. To understand this potential source of 
ambiguity in using the recursion relations at large N, consider the double-well potential, 
V{4>) = o~(j) + hfi 4> 2 + \g <^ 4 , with fx < 0, g > and a small symmetry breaking term a (p. 
To investigate symmetry breaking it is useful to study the effect of interchanging the limits 
iV — ► 00 and a — > 0± on the values of So = (O|0|O). By Z 2 symmetry, if we take a — >■ 0± 
followed by iV — > 00, we must get Sq = 0, whereas if we take iV — >• 00 followed by a — > 0± 
we have 

So = ±y^/g~, (2.3) 

as can be demonstrated by using steepest descent at the stationary minima of the poten- 
tial V(x). Similarly, to complete the necessary boundary conditions, one can show that 
Ri = (0\(j) 2 \0) — Sq takes on Ri = —fi/g and Ri = for these two limits respectively. 
Now we can in principle use the recursion relations to obtain both symmetric and bro- 
ken symmetry solutions. In general terms, this is just the familiar feature of spontaneous 
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symmtry breaking. For example the Ising model with no external field on a finite lattice 
must have (sj) = by symmetry, but in the large volume limit (at temperatures below 
the Curie point) the relevant (i.e. stable) solution is a broken vacua with non-zero values 
for (sj) = ±m obtained by applying an infinitesimal magnetic field. 

As we will see shortly, a more general possibility at infinite N is to characterize the 
"vacuum" state for the double- well <p 4 model by a mixing angle 9q, 

|0, O > = cos(6 /2) |0, +) + sm(0 o /2) |0, -) (2.4) 

where |0, ±) are orthonormal and 0|O, ±) ~ ±a/— /u/g|0, ±). As a consequence for the 
mixed state the boundary conditions becomes 

So = \/-v/9 cos(6 ), (2.5) 

instead of (|2.3p, with the constraint, 

R x = -fjt/g - S% > 0. (2.6) 



The goal of this paper is to understand how the solutions of the double-well matrix 
models depend on the initial boundary conditions, e.g., for the (p 4 model, on the first two 
moments, Sq and R±. In the large iV (or planar limit) we find a large class of solutions 
consistent with the broken symmetry boundary condition in addition to the (meta stable) 
pure gravity solution at higher free energy. Indeed we will show that the above qualitative 
discussion is a rigorous consequence of the planar solutions in the two-band ansatz. More 
generally however, we feel that the lack of a full understanding of the effects of this bound- 
ary condition represents an important gap in our ability to fully determine and solve the 
string equations resulting from matrix models. 
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3. Multiple Solutions In Matrix Models : Tree Level Analysis 

In this section we establish the existence of multiple solutions in two models: (i) the 
double well (f) 4 model and (ii) Gaussian Penner model. 



3.1. The Double Well <p A potential 

For the double well potential ( V(4>) = a<p + \ii 4> 2 + \g </> 4 ), Eq. (^2|) reduces to 
the recursion relations 

^ 2 2 

"TT = R n \P' + g(Rn+l + Rn + Rn-1 + S n + S n _ 1 + S n -lS n )\ 

N (3.1) 
= a + fiS n + g[R n+1 (S n+1 + 2S n ) + R n (2S n + 5 n _i) + S*]. 

We shall first illustrate our procedure by considering a symmetry-breaking solution under 
a period-one ansatz for both the R's and S's: R n — > R(jf), S n — > S(jf) ^ 0. Ignoring 
the 1/N corrections, Eq. ( |3.1|) leads to two relations, which allow us to solve for R(x) and 
S(x): 

R(x) = (1/15(7) [-» ~ Vv 2 ~ 15gx] 
S(x) = ±>J—n/g — 6R(x). 

With \i < and g > 0, one has -R(O) = 0, S(0) = ±y— Jijg, consistent with our discussion 
on symmetry breaking in Sec. 2, and R(x) monotonically increasing for < x < fi 2 /15g. 

The generating function F(z) = ^(^ r ~^) f° r a period-one ansatz at the tree level 
can be expressed in terms of R(x) and S(x) as follows: 

F(z) = dx — —— = / dx — (3.3) 

Jo Jo 2nz-<P(x,6) Jo ^[z-S(x)} 2 -4R(x) 

where 4>(x,9) = S(x) + a/ R(x)(e' ie + e~ lB ). One can verify by explicit calculation using 
(|3.2|) that the eigenvalues lie in an interval where z± = S(l) ± 2^R(1), with a 

single-band eigenvalue density, p(z), in agreement with the result of Shimamune (Ref. |l|) 
obtained by using the Schwinger-Dyson equation. On the line fx = —^/15g, R(x) develops 
a square-root type singularity at x = 1, leading to a Painleve-I equation in the double 
scaling limit, appropriate for the c = 2D gravity solution. 
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Since all the eigenvalues are concentrated in a single well, this symmetry-breaking 
solution does not correspond to a configuration with the lowest free energy, but is a sub- 
dominant solution. In the double scaling limit it is unstable against the tunnelling of a 
single eigenvalue into the other well, exactly like the subdominant solution for the pure 
gravity in the <f> 6 model. We now turn to solutions which have eigenvalues in both wells, 
which include a class of solutions for which p(z) is exactly Z2-symmetric (when a = 0) at 
the tree level, but which in general break this symmetry at higher order in 1/N. 



Let us consider a period-two ansatz for both the R's and S's 

71 71 

R n = A( — ), S n = C( — ) for n = even, 

7V 7V (3.4) 
R n = B(^), S n = D(^) for n = odd. 

Taking A, B, C and D to be continuous, and ignoring the 1/N corrections in the recursion 
equations (|3.1|) , we obtain four tree level recursion relations, 

2x = fJL ef f{A + B)+ g(A 2 + AAB + S 2 ), (3.5a) 

0= (A-BYA + B+ ^Z£), (3.56) 

9 

= 2a + n(C + D) + g[3(A + B)(C + D) + C 3 + D 3 } , (3.5c) 

= (C -D)(A + B+^), (3.5d) 

9 

where fj, ef f=fj, + g(C 2 + CD + D 2 ). 



Since we have already considered the single-band (pure gravity case) , we can assume 
that either A ^ B or C ^ D. To carry out an exhaustive analysis of the full set of 
solutions to these equations, it is useful to note that only three out of the four equations 
are independent. The three independent equations take the simple form, 

A + B -CD = - (C + D) 2 , AB = - , and V'[-(C + D)\ = 0. (3.6) 
9 9 

Note that there is no condition on C — D which can be independently chosen for every 
value of x. The first two equations allow one to express the explicit solution for A and B 
in the familiar form||, 

^=^(-^e//± v /^ // -4ox) B = -^(-(i eff =FyJij? eff -4gx). (3.7) 

In general for the double- well potential (for \i < and small enough a), there are three 
^-independent real solutions to V'[—(C + D)] = 0. With C + D fixed, it follows that the 
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first combination A + B — CD in (3J5) is also fixed. Consequently, a "circular" constraint 
oni-B and C — D can be found 

(A - B) 2 + — = 1[(C — .D) 2 + — + 3(C + -D) 2 ] 2 . (3.8) 

#16 # 



This constraint can be represented by contours for each fixed value of x, < x < 1, 
in a two-dimensional plane with A — Bos the vertical axis and C — D the horizonal axis. In 
fig. II we exhibit them for the class of solutions where a = and C + D = 0. The external 
contour, corresponding to x = 0, is precisely the constraint suggested in the qualitative 
discusion of symmtry breaking of Sec. 2. To see this one must consider carefully the proper 
definitions at the boundary, namely A — B = Rq — Ri = —R\ and C = —D = Sq at x = 0, 
which yields the contraint Eq. ( |2.6[ ) in the form, 

Rl = {Sl + ^/g) 2 . 

For x = x cr = /U 2 /4(7, the contour shrinks to a point, A — B = C — D = 0, about which a 



double scaling limit can be taken [TCj . 



The different solutions to (|3.6| ) can be parametrized by curves in this plane traversing 
from the x = to the x = 1 contours (see for example fig. 4(b)). For instance, a solution 
can be specified, at each value of x, by the polar coordinate, 6*(x), for the intersection of the 
curve with the contour (|3.8|) . Conversely, once 9(x) is chosen for every x, a unique solution 
to Eq. ( p.6| ), (A(x), B(x), C(x), _D(x)), is obtained |fT7| . In analogy with (|2.4| ), each state, 
\n), < n/N < 1, in the large N limit could be thought of as a linear conbination of 
"left-" and "right" states, specified by an arbitrary density function, 6{x). The choice of 
9(x) represents the ambiguity of solution at the tree level. (In fact the "orbit" need not 
even be continuous. However, if the orbit is discontinuous, derivatives of A, B, C, D are 
large and cannot be ignored, as is assumed in the tree level analysis.) 



An interesting feature of these broken symmetry solutions is that, within the class 
specified by one of the three values of C + D, they all have the same tree level eigenvalue 
density and free energy. To see this consider the generating function F(z) = jj(Tr ) 
for a general period- two ansatz at tree level |T8[| 



H*) ^ 1/2 /" ** Pf-^f' .... (3.9) 







vV - z(C + D) - (A + B - CD)] 2 - AAB ' 



Notice that Eq. ( |3.9|) involves precisely the three combinations, Eq. (|3.6|) , which are fixed 
by the four recursion relations. Therefore, once the value of C + D is chosen, one gets the 
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same generating function F(z). It follows that the tree level eigenvalue density p(x) and 
free energy are the same for all solutions within a class, i.e., are independent of the choice 
of9(x). 



For the most part we now shall restrict further discusions to the case a = and 
C + D = 0. All the solutions in this class give rise to the same eigenvalue density p(z). 
This class of solutions is continuously deformable into the limiting case A(x) 7^ B(x), 
C(x) = D(x) = 0, which is just the standard symmetric solution for the two-band solution, 
with p e ff = p. At the other extreme there is a maximally asymmetric two-band solution 
satisfying the condition A(x) = B(x), C(x) = —D(x) 7^ 0, with 

1/2 



A(x) = B(x) = R(x) = y f x/g, C{x) = -D(x)=± \p\/g-^Ax/g . (3.10) 

These values of C and D form the turning points at which our numerical solutions change 
the sign of A — B. (See fig. 1.) The symmetric solution corresponds to the choice 
6{x) = ±7r/2 and the maximally asymmetric solution to the choice 9(x) = 0, n. Our 
expanded class of solutions includes ones where the branch of the square root singularity 
(±) in Eq. ( |3.7| ) is exchanged between A and B as the trajectory rotates in the A — B 
versus C — D plane, as we note in the discussion of our numerical results for finite N (see 
Sec 3.4 and fig. 4 (c)). This corresponds to 9(x) winding around the circle a number of 
times as x goes from to 1. This is a precursor of the angular momentum variable of 
the double scaling solutions. The rigid implementation of Z<i symmetry and the boundary 
conditions on the recursion coefficients would have yielded only the 9(x) = ±n/2 solution. 
All other solutions correspond to a breaking of the Z<i symmetry. 



3.2. Gaussian Penner Model 



The second example, we would like to consider, is the Gaussian Penner model. The 
potential for a general Penner model is V((p) = Vo(</>) —t log <fi, where Vq is a polyno- 
mial[jl2|. If Vq{4>) = 4> the model is the linear Penner model |l9(|p0| , if V (4>) = p 4> 2 /2 the 



model is the Gaussian Penner model |TI| [|13|, where we interpret the log term as | log </> 2 . 
In fig. |^ we display the Gaussian Penner potential for different values of t. Consider first 
the situation where t > 0. (The region t < is reached by analytic continuation. p4| ||12|| ) 
As the potential is a double well, the period-two ansatz may be applied here also. 



The recursion relations (Eq. ( |2.2|) ) for a general Penner model reduce to 



fry , , 

— = v^(n-l|W)|n)- V^(n-l|0~V> 
= (n\Vo((f>)\n) -t(n|^ _1 |n)- 



(3.11a) 
(3.116) 



Let us denote W n = y/R n (n — l|V^(0)|n) and Y n = (n\VQ(cf))\n) for later notational conve- 
nience. For the Gaussian Penner model, W n = fxR n and Y n = fiS n . 



Eqs. ( |3.11a, fc| ) are unusual since they involve matrix elements of 4> . For t > 0, 
they can be solved in the spherical limit by a procedure similar to that used for deriving 
the generating function F(z), Eq. (fT9l), under a period-two ansatz plR . By considering 
n even and odd, Eqs. ( |3.11a, b\j should normally lead to four conditions. Just like the </> 4 
model, only three are independent, and they can be cast in the following form: 



C + D = 0, 

2x + t 



A + B-CD 



AB 



x(x + 1) 



(3.12a) 
(3.126) 

(3.12c) 



For the symmetric solution where C = D = 0, one finds 



A(x) 



x 



B(x) 



x + 1 



For the maximally asymmetric solution, on the other hand, one has 



(3.13) 



A(x) = B(x) = -Jx(x + t), C(x) 2 = - 
fj, n 



(2x + t) -2y / x(x + t) 



(3.14) 



Observe that Eqs. (|3.12a — q) are precisely the necessary combinations which are 
needed in Eq. ( |3.9|) for determining the generating function of our symmetric Gaussian 
Penner model in the spherical limit leading to symmetric two-band structure. Therefore, 
for this class of solutions and in particular for the symmetric and maximally asymmetric 



solutions, the eigenvalue density and free energy are again identical at tree level [22 



3.3. Correlation Functions that distinguish between symmetric and asymmetric ansatz 
solutions 



Given that the tree level generating function and hence eigenvalue density and free 
energy are the same, one might ask if there are other correlation functions that distinguish 
between the various solutions. It is normally assumed that after taking the period-two 
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ansatz the large N limit is smooth for all correlators. However, consider the correlator 
(Tr<f)Tr<fi) c . In terms of recursion coefficients, 



(Tr^Trcj)). = R N . (3.15) 

In the symmetric solution, since the Rn alternate between An and Bn as N goes from 
odd to even, this correlator at large N depends on whether oo is approached through odd 
or even N. In one case it is -4(1), in the other B(l). For the </> 4 model, they differ even 
at tree level by a//U 2 — Ag/g which is singular at criticality, (see Eq. (|3.7|) ). On the other 
hand in the maximally asymmetric solution Sn is period two but Rn is of period one, 
(see below Eq. (|3.10|) ) hence this correlator has no discontinuity between odd and even. 
This is an example of a correlator that distinguishes between the two solutions. A similar 
difference between odd and even N is known to exist in the context of unitary matrix 



models [23 . 



Another example is (Tr<pTr<pTr(j)) c = Rn (Sn-i — Sn) • In the symmetric solution 
Sn = and this vanishes. But in the maximally asymmetric solution since S n is period 
two, (and C = —D), Sn-i — Sn changes sign as one goes from odd to even N. In particular, 
(Tr<f)Tr<j)Tr<j)} c = (\n\/g 2 — 2(l/(?) 3 / 2 ) 1 / 2 for N odd and its negative for N even. 

Another characterization of the difference between the solutions is the following: If 
one truncates the infinite dimensional matrix (n\(p\m) to an N x N matrix corresponding 
to the subspace of the first N orthogonal polynomials, the eigenvalues of this N x N 
matrix are a good approximation to the saddle point configuration of the eigenvalues at 
large N. Since its matrix elements are given in terms of R's and S's, we can determine 
the eigenvalues numerically from a given solution of the recursion coefficients. In fig. [|, 
we show the locations of the eigenvalues so obtained for the tree level symmetric and 
maximally asymmetric solutions for N = 24, 25. For N = 24, half the eigenvalues are in 
one well and half in the other, for both solutions. However, for N = 25, (odd A), there is 
a striking difference between the two solutions. For the asymmetric solution, there is one 
extra eigenvalue located in one of the two wells, (the well selected depends upon the sign 
of So ) , whereas for the symmetric solution, this extra eigenvalue in the the center (on top 
of the barrier), thus preserving the symmetry between both the wells. 



3.4- Numerical Approach to Finite N Solutions 



Another approach to understanding the role of these multiple solutions is to pursue 
a numerical study of finite N solutions and attempt to take N large enough to see a cross 
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over to the large iV (or double scaling) regimes. Although we will postpone a detailed 
analysis of our results, there are several general features which can help to understand the 
present discussion. The recursion relations for the double well potential follow from the 



variation of an effective action, [124" 



V ef f(Rn, S n ) = + (iRn + | (i£ + 2R n R n+1 ) 

n=0 

i C , ^C2 , 9q4 ( 3 - 16 ) 
+ gRn{Sn + S^-i + S n -iS n )}, 

with the defintions, S-i = and Rq = 0. Morover we must take g > 0, if the effective 
action (like the actual potential) is to be bounded from below. This formalism provides 
a natural way to investigate our set of multiple solutions, by removing the boundary 
conditions on £0 and R\ and replace them with the asymptotic condition that R n and 
S n are smooth functions as n — > 00. This is simply the one band ansatz in the extreme 
limit of n/N very large. Therefore it is again interesting to ask what is the full set of local 
minima. 



Earlier work on symmmetric solutions for the degenerate three-well potential have 
observed the recursion coefficients with very complicated, "chaotic looking" behavior, when 
calculated by a numerically method logically equivalent to minimizing an effective poten- 
tial |25|. We also have observed complicated behavior for two degenerate wells when we 



allow symmetry breaking terms (S n 7^ 0), which we have been able to relate to the exis- 
tence of our multiple solutions in the planar limit. As an illustration consider the solution 
presented in fig. f|. However due to the degeneracy of multiple solutions at N^, great care 
must be taken with the minization procedures. 



For example we have minimized V e ff for the double- well potential with N = 512, 
fj, = —2, g = 1 and a = 0.1, starting from a random distribution of 2048 coefficients for 
R n and S n . Using a variety of minimizaton procedures on the CM-5 at Boston University 
and the NeXT station at CERN, we see that after only several 100 iterations the curves 
conform roughly with the large iV constraints but they can have a great variety of coutours 
in the A — B vs C — D plane. However if we go further for another 100,000 iterations, a 
smooth spiral curve (see fig. 4 (c)) begins forming near critical x [x cr = 1), with the large 
iV constraints improving to about 1% as might be expected in a transition region from 
one of our N = 00 solution to a particular double scaling solution with non-zero orbital 
quantum number / 7^ 0. After 100,000 iterations the value of V e ff/N departs from its 
theoretical N = 00 value by 0.0034. The final results on questions as to the stability of 
non-zero orbital solutions, the possibility of residual degeneracies at finite N and especially 
the existence of choatic regimes require accurate and non-trivial compuatational power. 
Further details on this as well as a study of higher 1/N corrections will be presented in a 
future publicationp6[]. 
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4. The Double Scaling Limit 



4-1. Double Well 4> 4 Potential 



The double scaling equations for the <p A model have been discussed by a number of 
authors PJ2^|5|-|IU| . The steps involved in the double scaling analysis of symmetric breaking 



solutions for a Z2 symmetric potential are the same as that for solutions of a general 
asymmetric potential |]7|-|T0[| , since in both cases one includes both R n and S n in the analysis. 
One sets x = 1 — t 2 t, (recall x = n/N) and e = iV -1 / 3 . For the symmetric solution, 
C n = D n = 0, while A n and B n are 

A n = a + e(/ e (t) + f (t)) + e 2 (r e (t) + r (t)) + (4.1a) 

B n = a + e(f e (t) - f (t)) + e 2 (r e (t) - r (t)) + (4.16) 



On substituting this symmetric double scaling ansatz into the recursion relations (| 
and equating terms with powers e°, e 1 , e 2 , ande 3 , we get eight equations, two of these are 
used up by a Q (e° equations), the tree level result. (Note a a = —fi/(2g); in what follows, 
we adopt the convention where \x = —2 and g = 1). That leaves us with six unknowns 
and six equations hence all the unknowns can be determined. Most of them are zero 
(e.g. f e = r a = 0....), while the others are determined in terms of f (t) = f(t), e.g., 
r e = (f 2 — £)/4. The function f(t) satisfies the Painleve-II equation 

f" ~ \ f + I ft = 0- (4-2) 
The suseptibility X ~ ~ / 2 /2 - r e = (/ 2 + t)/4. 

For the maximally asymmetric solution, the double scaling ansatz for C n and D n 

are 

C n = eg(t)+e 2 (4.3a) 

D n = -eg(t)+e 2 (4.36) 

and A n = ao + e 2 r e (t) + , B n = a D + e 2 r e (t) + • • •. Substituting this maximally 

asymmetric double scaling ansatz into the recursion relations and equating powers of e we 
get, r e = — (g 2 + £)/4, etc., and 

9"-\g 3 + \gt = 0, (4.4) 
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the same as Eq. (|4.2|), with g replacing /. Under this ansatz, the suseptibility x can be 
expressed as (g 2 + t)/4. 



We next consider the general symmetry breaking solutions, where A — B ^ and 
C — D 7^ in the planar limit. Substituting the double scaling ansatz, Eqs. (|4.1| ) and ( fO| ) 
into the recursion relations and equating powers of e we get r e = (f 2 — g 2 — £)/4, etc., and 
the following coupled equations [f7|-|10| 

f"-f(g 2 + f 2 )/4 + ft/2 = 0, (4.5a) 
g" -g(g 2 + f 2 )/4 + gt/2 = 0. (4.56) 

The suseptibility is now given by % = (f 2 + 9 2 + /4- 

While Eqs. ([4.5|) have been obtained previously in the context of asymmetric poten- 
tials, we would like to emphasize that they describe the multiple (and in general symmetry 
breaking) solutions that exist even for a ^-symmetric potential. To see the symmetry 
breaking nature of these solutions more explicitly and to make contact with the tree level 
discussion in the previous section, introduce a two-dimensional vector r = (g, /), in terms 
of which the coupled equations can be written as r — (l/4)(r 2 — 2t)f = 0. We can next 
make a change of coordinates 

f = r sin8(t), g = r cos8(t), (4-6) 

so that x ~ (r 2 + 1)/4 and the coupled equations become 

11/2 

f r 3 + - rt o=0, (4.7a) 

4 2 r 6 

r 2 8 = I. (4.7b) 
Note that since A — B oc / and C — D oc g in the double scaling limit, the variable 8 is 



the same as that introduced in the previous section. The constant / in Eq. (|4.75|) is the 
"angular momentum", and it is a constant of the motion due to the U(l) invariance of 
Eqs. (|4.5|) . Note that for / = equation ( |4.7a|) is just the Painleve-II equation in the r 



coordinate. Thus in the / = sector of this model, we reproduce the same double scaling 
results for both the symmetric and maximally asymmetric solutions. But for the / ^ 
sector the double scaling equation is different; hence the behavior of the system in this 



sector for the multiple solutions is different. Iterating Eq. ( 4.7a ) at large t, one finds 



x =|-(i±i!)r 2 + ... (u 
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At the tree level we had a degeneracy of solutions parametrized by the function 0(x), 
whose value could be independently chosen for xe[0, 1]. The double scaling analysis based 
on the (4.1 )and (4.3 )tells us that only a two parameter family of these solutions, labelled 
by I and one global rotation angle 9q, survives in the double scaling limit. (6(x) is no longer 
any function of x, but constrained such that r 2 9 = / is a constant.) Solutions labelled by 
different values of I give rise to the same susceptibility at tree level (the first term in ( |4.8| ) 
is /-independent) as expected, but differ at higher orders. 



A relevent question that arises is: Just as the symmetric solution (S n = 0) is the 
"natural" solution that follows from the symmetric potential V(x), ("natural" in that it 
respects the Z<i symmetry of V(x), and the recursion coefficients are specified by initial 
conditions given by the integrals discussed in section 2), is there a perturbed potential 
of which the symmetry breaking solutions are natural solutions? This is presently being 
investigated. At this point we remark that the perturbation cannot be a rigid translation 
of the potential, which induces linear and cubic terms. Such a perturbation has been 
discussed in Ref. It can be easily seen that the natural solution to the shifted 

potential Vb(x) = V(x — b) is given by the same R n as for V(x), with S n = b = const. 



Should one decide to introduce a small explicit symmetry-breaking term, a = 
N~ 2 / 3 a, as was done in Ref. 0, the vector equation above remains U(l) invariant by 
simply adding a "constant magnetic field" term of the form B^r x z, where Bq oc a and 
z denotes a unit normal to this two-dimensional plane. One can thus again reduce it to a 
single radial equation with the solution depending on a "generalized angular mumentum" , 
I = r 2 9 + B r 2 /2 = constant, 



B Q r 2 \2 



1 o 1 (I- ., , 

r- - r 3 + - rt- - ^—^ = 0. 4.9 

4 2 r 3 



4-2. The Gaussian Penner Model 



The double scaling solutions for the Gaussian Penner model have been discussed in 



Refs. |T2[ and [T^]. We reproduce the proofs below for completeness (also simplifying them 
somewhat). The critical point is t = —1. By strictly enforcing symmetry, this model 
can be solved exactly first at t > 0, so that the criticality at t = —1 can be exhibited 
explicitly. Note that (?t,|(/> _1 \n) = by Z 2 symmetry, Eq. (|3.116| ) thus reduces to S n = 0, 
the symmetric ansatz. Eq. fl3.11q|) can also be solved exactly for all n. Since W n = fxR n 
and since In — l|^ _1 |n) = for n even, (n — l|^ _1 |n) = -7= for n odd, it follows that, for 
n even, R n = n/fiN, and for n odd, R n = (n + t/V")///iV"||12||. 
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Since we know the exact result for the R' n s for the symmetric solution, the exact free 
energy may be obtained 



N/2-1 

r= k log[(2k + fx+l) (2k + n-l) 

k=i 



(4.10) 



where t = —1 + On expanding the free energy in powers of \x we get 



r = - ii 2 log fi + ^_ i g n _|_ .... 



(4-11) 



The coefficient of the second log/U, xi — 1/12, comes from the torus contribution, indicates 
that this free energy cannot be identified with the Legendre transform of the free energy 



of the c=l string at self dual or twice the self dual radius [12 



It has been stressed in Ref. [|12] that the exact solution to the Gaussian Penner 
model is characterized by the fact that B(x) has a linear zero at x = 1 when t = — 1 
while A(x) is non-zero there, (see Eq. ( |3.13|) ). This same feature also holds in general 
for symmetry breaking solutions where C(x) = —D(x) ^ 0. However, the maximally 
asymmetric solution provides an exception to this rule. When A{x) = B(x), one has 
A(x) = B(x) ~ (x + t) 1 / 2 , (see Eq. (|3.14|) ). That is, both A(x) and B(x) have square-root 
type behavior near x = 1 in the spherical limit at t = —1. Since it is the behavior of 
R n near x = 1 which determines the criticality of the model, it follows that the resulting 
double scaling limit for the maximally asymmetric ansatz could be non- generic. 



Let us next concentrate on the maximally asymmetric solution. We note first that, 
with S n ^ 0, Eqs. ([ma|) and ( $7\TQ ) can be re-written as |T| [0 



W n + W n+1 + S n Y n = 2U + 1 N +N \ (4.12a) 
S N [W n+1 -W n -^)= R n Y n _ x - R n+1 Y n+u (4.126) 

so that matrix elements of <p would not appear explicitly. For the maximally asymmetric 
solution the double scaling solutions may be found as follows. With x = 1 — z/N, t = 
— 1 + /jl/N, e = 1/x/iV, A n , C n , and D n can be expanded as []13| 



A n = l/VN p(z) + ... 

C n = l + 1/VN a(z) + 1/N (7i (z) + .... (4.13) 
D n = - (l + 1/y/N a{z) + 1/N a x {z) + ....) . 
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On substituting into the recursion relations (4.12 )and equating equal powers of N, 
we can determine cr(z) in terms of p(z) and o~i(z) in terms of o~(z). The equa- 
tion for p is p(z)p(z — 1) = p — 1/2 + z and the double scaling result for R n is 
R n ~ T (§(iV - n + // + 3/2)) /r (|(iV - n + p + 1/2)). The double scaled free energy 
is T = J2k=i 1 ^ /o^ [(2/c + ^ + l/2)(2/c + — 1/2)] plus/i independent terms. On ex- 
panding in powers of /x one gets 

1 5 
r = - p 2 log p - — log p (4.14) 

The coefficient xi = —5/48 confirms that this criticality corresponds to that for the free 
energy of the c=l string at twice the self dual radius as conjectured in |2lf . We note that 



although the tree level free energy gave identical results for the symmetric and maximally 
asymmetric ansatz (see remarks below Eq. ( |3.14j )), the double scaled free energies ( f4.11| ) 
and ( f4.14j ) are quite different. 
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5. Discussion 

The existence of multiple double scaling solutions at the same critical point, which 
share the same tree level behavior, is an unusual and previously unnoticed feature. We 
have exhibited this behaviour in two completely different models-the double- well <p 4 po- 
tential and the Gaussian Penner model. Both these models possess Z2 symmetry, and a 
charateristic feature of the class of solutions we consider is the fact that only one of the 
solutions respects this symmetry to all orders, all other solutions break this symmetry. 
Although the physical consequences of symmetry breaking in matrix models are not yet 
fully understood, it is useful to consider the analogy of the multiplicity of solutions here 
to the property of coexisting ferromagnetic phases below the Curie point. In the Ising 
model there is an infinite set of mixed phases with identical free energy per unit volume 
in the infinite volume limit. The domain walls that characterize the mixed phases give 
rise to lower order contributions in the free energy expanded around the infinite volume 
limit. The "bulk" contribution in the multiple solutions is the same (the tree level eigen- 
value density is the same) but they differ by amounts suppressed by powers of 1/N, like 
"surface" contributions. 

A significant difference between the Gaussian Penner model and the symmetric </> 4 
model is that in the Gaussian Penner model there does not seem to be an angular momen- 
tum parameter / characterizing the double scaling solutions. Further, the double scaled 
free energies of the symmetric and maximally asymmetric solutions already differ, unlike 
in the </> 4 case where both of these were / = solutions with the same double scaled free 
energy. In spite of these differences, however, both models display the same general phe- 
nomenon, namely, that the enlarged class of symmetry breaking solutions produces the 
same free energy at tree level, and contain solutions that produce different free energies at 
higher orders. We expect this to be a generic feature of multi-cut matrix models. The exis- 
tence of multiple solutions is related to the fact that when the potential has more than one 
minima, the number of smooth functions required to represent the recursion coefficients 
exceeds the number of constraints obtained from the recursion relations. Thus multiple 
solutions will exist even when the potential has no symmetry as we noted for a ^ in the 
</> 4 model. 

It would be interesting to know whether and how these solutions can tunnel into 
each other. This would be particularly interesting for the Gaussian Penner model where 
the maximally asymmetric solution corresponds to the c = 1 string compactified at twice 
the self-dual radius. We are now studying in greater detail higher order terms beyond the 
planar approximations as well as numerical solutions at finite N to the effective potential, 
Eq. ( pi.l6| ), to determine more precisely the relationship between multiple solutions of the 
planar versus the double scaling limit. 
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Figure Captions 



Fig. 1. Constraint on A — B versus C — D for x = 0, 0.5, .75, .875, .95, with \i = — 2 and 
9 = 1- 



Fig. 2. The Gaussian Penner potential for t < 0, t = and £ > 0. 



Fig. 3. Eigenvalue distributions for even and odd iV: (a) iV = 24, symmetric solution; 

(b) iV = 24, asymmetric solution; (c) iV = 25, symmetric solution; (d) iV = 25, 
asymmetric solution. 



Fig. 4. Graphs of recursion coefficients for a sponteneously broken solution of the double- 
well potential, (a) The R n and (b) the S n coefficients after 100,000 minimizaton 
steps from a random start, (c) Orbit in the A — B vs C — D plane. 
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